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The upcoming reactor neutrino experiment, KamLAND, has the ability to explore the Large 
Mixing Angle (LMA) solution to the solar neutrino problem. Here, we investigate the precision 
to which KamLAND should be able to measure these parameters, utilizing the distortion of the 
energy spectrum of reactor neutrinos. Incomplete knowledge of the fuel composition of the reactors 
will lead to some error on this measurement. We estimate the size of this effect. 
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I. INTRODUCTION 

For a long time we have been aware that there is a 
deficit of neutrinos emanating from the sun. An apparent 
solution to this problem is the phenomenon of neutrino 
oscillations Q. Particularly in light of the compelling 
atmospheric neutrino data from SuperKamiokande 0, 
we fully expect solar neutrinos to oscillate. Oscillations 
between two flavors can be described effectively by two 
parameters: the mass difference, Am 2 , and a mixing 
angle, tan 2 9. There are four regions in the Am 2 , tan 2 
plane that explain the observed data on neutrinos from 
the sun. KamLAND is designed to explore one of 
these regions, the Large Mixing Angle (LMA) solution, 
through the detection of reactor neutrinos. 

KamLAND |j| is an experiment to be located at the 
old Kamiokande site in the Kamioka mine in Japan. This 
location is of key importance to this experiment, as it is 
situated in the vicinity of 16 nuclear power plants which 
will contribute a significant neutrino flux. KamLAND 
consists of approximately 1 kiloton of liquid scintillator 
that will detect reactor neutrinos through the reaction: 



P- 



(1) 



The positron is then detected when it scintillates and 
when it annihilates an electron. This annihilation, in 
delayed coincidence with the 7-ray from neutron capture, 
represents an easily recognizable signal. 

In this letter, we assume that the solution to the solar 
neutrino problem is in the LMA region. We explore 
the accuracy to which the KamLAND experiment can 
utilize the reaction of Eq. (|l|) to measure the param- 
eters of a LMA solution. In section [n] we review the 
basic procedure for computing the spectrum of neutrinos 
created at reactors. In section III, we review how to 
compute the expected number of events at KamLAND. 
In IV, we describe the results of an analysis of the energy 
spectrum of the detected neutrinos. Next, in section |v| 
we address the question of systematic errors associated 
with an incomplete knowledge of the incident neutrino 
spectrum. Finally, we conclude with a brief discussion 
of other possible systematics to be investigated. We 
also briefly mention some of the implications of this 
measurement for future neutrino experiments. 



II. DETERMINATION OF ENERGY SPECTRUM 

Because the mixing angle of v e in the mass gap 
responsible for the atmospheric neutrino oscillations is 
constrained to be small from CHOOZ reactor neutrino 
experiment [Q, the solar and atmospheric oscillations 
decouple to a very good approximation. In this approxi- 
mation, the important equation for the analysis of reactor 
neutrino oscillation becomes: 



sin 2 20 sin 2 



1.27Am 2 (eV 2 )L(km)* 
. E(GeV) - 



(2) 



One of the advantages of the KamLAND design is that 
it is expected to have good energy resolution. An energy 
resolution of ^jp- = or better is anticipated, where 
E is measured in MeV ||. As a result, if there are 
oscillations, one might hope to utilize the E dependence 
of Eq. (^) to assist in making an accurate measurement 
of the oscillation parameters. 

Of course, in order to take advantage of this energy 
dependence, one must have knowledge of the energy 
dependence of the incident neutrinos. That is to say, one 
is required to have a good knowledge of the unoscillated 
spectrum. Wc discuss the determination of this spectrum 
here. 

A number of short baseline experiments || have mea- 
sured the energy spectrum of reactors at distances where 
oscillatory effects should be completely negligible. A 
phenomenological parameterization of these spectra ex- 
ists H, which depends on the isotope involved. Namely, 
Vogel and Engcl find that: 



dN v 
dE v 



c ,a Q +a 1 E„+a2E v 



(3) 



where the fitted values of the parameters at are repro- 
duced from H in Table |j| This spectrum is given in 
units of „,„"' c „:- ■ Therefore, given this spectrum, it 



McV— fission * 

remains to determine how many fissions of each isotope 
there are. For a given reactor site, this will depend on 
three factors: 

1. The thermal power of that reactor. 
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TABLE I. Parameters for ^f- parameterization, 
resulting spectrum is given in units of 77 e /MeV-fission. 



The 



Isotope 




^ 9 Pu 


238 yj 


Pu 


ao 


0.870 


0.896 


0.976 


0.793 


ffli 


-0.160 


-0.239 


-0.162 


-0.080 


(12 


-0.0910 


-0.0981 


-0.0790 


-0.1085 



2. The isotopic composition of the reactor fuel. 

3. The amount of thermal power emitted during the 
fissioning of a nucleus of a given isotope. 
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FIG. 1. The time dependent composition of reactor fuel. 
Taken from reference ]9|. 

For convenience, we reproduce the relevant reactor 
characteristics in Table |TJ. The maximum thermal power 
is given in the table. This addresses the first point above. 
Isotopic composition is a somewhat complicated issue, 
and we will investigate it more fully in section [v| For 
the present, let us simply note that, in general, the time 
composition of the fuel varies roughly as in Fig. [l]. For 
the analysis of section IV, we take all reactors as having 
a composition varying as in this figure. We assume that 
at the end of the 7 month cycle shown, the composition 
reverts back initial levels (refueling). Finally, for point 
three, we note that the thermal energy, associated 
with the fissioning of a given nucleon is known [j7|. 



The relevant values are displayed in Table III. With 
these pieces of information, we are able to determine the 
neutrino spectrum emitted from each reactor, S(i,t,E). 
Here, i labels the reactor. 



III. DETERMINATION OF EXPECTED 
NUMBER OF EVENTS 

Now that we have the initial spectrum in hand, we 
review how one finds the expected number of events 
at KamLAND. To determine the number of neutrinos 
detected at KamLAND, one must convolve the cross 



TABLE II. Reactor parameters. Reproduced from refer- 
ence tel. 
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TABLE III. Energy per fission of isotopes that make 
significant contributions to the thermal power of a reactor. 



Fissioning Isotope 




^ y Pu 


•2 l M jj 


M1 Pu 


Energy Per Isotope (MeV) 


201.7 


205.0 


210.0 


212.4 



section, a(E v ), for the reaction shown in Eq. (Q) with 
the reactor spectra, S(i,t,E). To lowest order the cross 
section is 



o{E v ) 



2tt 2 



m\fr n 



PeE e . 



(4) 



Here / = 1.69 is the integrated Fermi function for 
neutron /3-decay, m e is the mass of the electron, E e is 
the electron energy, and p e is the electron momentum. 
The energy of the electron, to lowest order, is given by 



E e =E U - 1.293 MeV. 



(5) 



That is to say, the energy of the electron is basically the 
energy of the incident neutrino minus the proton-neutron 
mass difference. In our numerical calculations, we used 
the cross section which takes into account the nucleon 
recoil, which may be found in pd| ]. 

The expected number of events at KamLAND is given 
by: 



N(t,E„, Am 2 , sin 2 29) 



E 



S(i,t,E u ) 



a (E V )P (Am 2 , sin 2 20). 



(6) 



Here, P is the probability from Eq. (g), u(E v ) is the 
full cross section analogous to Eq. (|), and S(i,t,E y ) is 
the initial energy spectrum, to be calculated as described 
above. The values of the distances can also be found in 
Table || 
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IV. RESULTS 

With the expected number of events in hand, it is only 
a matter of a simple \ 2 analysis to fit for the oscillation 
parameters. In this section, we ignore background effects. 
To estimate the precision to which KamLAND could 
measure the oscillation parameters, we assumed a 3 kt- 
year exposure where all reactors operated at 78% of their 
maximum capacity. It was assumed that the fuel compo- 
sition of each reactor varied as shown in Fig. |l| Perfect 
detector efficiency was assumed. To lowest order, this 
is not a bad assumption, given the recognizable delayed 
coincidence signal. The measurements that KamLAND 
could be expected to perform with these assumptions are 
shown in Fig. g. Contours were generated by finding the 
minimum % 2 , and then calculating the confidence levels 
for two degrees of freedom (tan 2 9 and Am 2 ). Data were 
binned in 0.5 MeV bins. The fit was done for visible 
energies^ between E mslUe = 1.22 MeV and E msMe = 7.0 
MeV. The first bin is smaller, since our parameterization 
of the incident spectra is only good above E v = 2 MeV. 
The cutoff at high energies is to avoid the low statistics 
bins, where there is not much statistical discrimination, 
and Poisson statistics would be necessary. 

Since the oscillations in this case depend on the mixing 
angle only through sin 2 29, there is a two-fold degeneracy 
in the measurement (hence the reflection symmetry 
about tan 2 9 — I). However, the LMA solution, which 
is overlayed, does not posses such a symmetry, so it is 
necessary to plot against tan 2 9 and not sin 2 29 12 1. 

From the figure, it is clear that KamLAND is able 
to make a very accurate measurement of the Am 2 
parameter in particular. To see why this is so, it is 
instructive to look at the quantity: 



10" 



R 



#Observed 
#No Oscillation 



(7) 



Plots of this quantity are shown in Fig. || This should 
essentially be the oscillation probability. By determining 
the position of the dip in the oscillation, one is able to 
determine the value of Am 2 . From the figure, one can 
see that the location of the dip can be well determined by 
KamLAND if the oscillations parameters lie in the LMA 
solution to the solar neutrino problem. 



V. FUEL COMPOSITION EFFECTS 

Let us now address the question of the effect of the fuel 
composition. To calculate the evolution of the different 
components of reactor fuel is a complicated business. 
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*Here, visible energy is defined as the energy seen in the 
detector from the e + and its annihilation, E visi i,i e — E e + + 
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FIG. 2. The expected measurement of Am 2 and tan 2 6. 
Contours for 68%, 90% and 99% CL are shown.The LMA 
solution to the solar neutrino problem is overlayed [^j. 

KamLAND expects to obtain the time dependent compo- 
sitions from the power companies involved. This section 
serves to demonstrate that this is an important piece of 
information for this experiment. 

To estimate the effect that an incomplete knowledge 
of the fuel composition would have on the measurement, 
we took the extreme cases. We look at the neutrino 
spectrum that would result from taking the composition 
to be as illustrated in Fig. |2|, we refer to this as the 
"true composition." We then tried to fit the observed 
data that were generated from the true composition 
with the expected number of events from three different 
compositions: 

1. The true composition, 

2. A composition that stays constant and equal to the 
composition at t = 0, 

3. A composition that stays constant and equal to the 
composition at t = 220 days. 

The different spectra are shown in Fig. ||. Although the 
basic shapes are relatively consistent, there are important 
differences, as can be seen by examining Fig. where we 
plot the difference between the spectra, normalized to the 
average spectrum. At the higher energies, differences in 
the initial and final spectra can reach over 15%. 

The difference between the three fits done in this 
way should give an estimate of the systematic error 
involved. A contour plot showing the measurements in 
the (tan 2 9, Am 2 ) plane is shown in Fig. ||. Although 
the values of the measured quantities are not drastically 
affected, particularly in the cases of moderate Am 2 , fuel 
composition clearly is an important systematic error to 
be considered. It is also of interest to see how incorrect 
assumptions about fuel composition degrades the fit. We 
show the x 2 f° r 10 degrees of freedom in Table IV. In 
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FIG. 3. The ratio, R, of the number of events seen, to the 
number of events expected if there were no oscillations. The 
error bars are statistical. 




E [MeV] 

FIG. 4. Neutrino spectra generate using the initial, final 
and average fuel composition. Although the differences here 
are somewhat difficult to discern, they are move easily visible 
in Fig. B. 
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FIG. 5. Here we have plotted the percentage difference in 
the initial and final spectra. More precisely, we have plotted 
the 100 times the difference between the initial and final 
spectra, divided by the average spectrum. 
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particular, the x 2 is not nearly as good for the fits with 
the final spectrum. 
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FIG. 6. 90 % confidence level contours for measurements of 
Am 2 and tan 2 8. The contours were generated by varying the 
assumptions about the composition of the fuel. For further 
explanation, see the text. 



VI. FLUX UNCERTAINTY EFFECTS 



TABLE IV. The fits corresponding to three different fuel 
compositions. For each point in (Am 2 , tan 2 (9) space, the 
events incident on KamLAND were generated with a time 
dependent (T) composition. Fits were done with the time 
dependent (T), initial (I), and final (F), composition. 



Comp. 


A 9 

Am 2 


9 ^ 

tan 2 9 


Fit Am 2 


Fit tan 2 


X 2 (10 dof) 


T 


1 x 10~ 5 


0.52, 1.92 


9.5 x 10 _c 


0.64, 1.55 


12.89 


I 


1 x 10~ 5 


0.52, 1.92 


1.1 x 10~ b 


0.43, 2.30 


9.64 


p 


1 x 10 -5 


52 1 92 


8.7 x 10 _o 


80 1 24 


21.67 


T 


2 x 10~ 5 


0.52, 1.92 


1.9 x 10~ b 


0.54, 1.87 


12.42 


I 


2 x 10~ 5 


0.52, 1.92 


2.0 x 10~ b 


0.53, 1.87 


9.50 


F 


2 x 10~ 5 


0.52, 1.92 


1.9 x W' b 


0.50, 2.01 


18.26 


T 


3 x 10" 5 


0.52, 1.92 


3.0 x W' b 


0.50, 2.01 


12.17 


I 


3 x 10~ a 


0.52, 1.92 


3.0 x 10~ ; ' 


0.51, 1.94 


10.2 


F 


3 x 10~ 5 


0.52, 1.92 


3.0 x 10~ b 


0.46, 2.16 


16.19 


T 


7 x 10" s 


0.52, 1.92 


7.2 x 10~ a 


0.45, 2.23 


10.73 


I 


7 x 10~ s 


0.52, 1.92 


7.0 x 10 _t> 


0.50, 2.01 


11.96 


F 


7 x 10~ 5 


0.52, 1.92 


7.2 x 10~ b 


0.39, 2.58 


11.57 


T 


1 x 10-* 


0.52, 1.92 


1.0 x 10~ 4 


0.51, 1.94 


12.96 


I 


1 x 10 -4 


0.52, 1.92 


1.0 x 10 -4 


0.58, 1.72 


8.58 


F 


1 x 10 -4 


0.52, 1.92 


1.0 x 10 -4 


0.45, 2.23 


21.87 



Finally, we relax our assumption that KamLAND will 
have complete knowledge of the overall flux normaliza- 
tion. We assume that the flux is known with a three 
percent error, consistent with the size of the errors in the 
/3-spectroscopy experiment at the Gosgen reactor (9| . Wc 
then can compute a x 2 , that takes this uncertainty into 
account. In particular, we write: 

2 -fa (N-l.0) 2 \ 
X = mm Xuns + 2 • ( 8 ) 

V U N J 

Here, a% — .03 2 . We vary the possible normalizations, 
N, and take the minimum of the right hand side. Of 
course, the xlins ^ s a function of N as well, as the 
theoretical prediction depends on the incident flux. Note, 
in general, this procedure will provide a smaller % 2 , at 
each set of oscillation parameters, resulting in larger 
contours. From Fig. [?], we can see that the contours 
are enlarged slightly. In particular, the uncertainty 
in the overall flux causes a slight degradation of the 
measurement of the mixing angle. This is a well-known 
and important effect. For completeness, we tabulate the 
values of \ 2 m Table [v| 

VII. EXPERIMENTAL SYSTEM ATICS 

Although we are not in a position to make a quantita- 
tive study of the experimental systematics at KamLAND, 



TABLE V. The fits corresponding to three different fuel 
compositions. For each point in (Am 2 ,tan 2 #) space, the 
events incident on KamLAND were generated with a time 
dependent composition. Fits were done either allowing the 
overall flux to float (Float) , or keeping it fixed at the predicted 
value. (Fixed). 



Comp. 


Am 2 


tan 2 9 


Fit Am 2 


Fit tan 2 9 


X 2 (10 dof) 


Fix 


1 x 10 _s 


0.52, 1.92 


9.5 x 10"" 5 


0.56, 1.55 


12.89 


Float 


1 x 10 _b 


0.52, 1.92 


1.0 x 10" 5 


0.56, 1.79 


12.38 


Fix 


2 x 10 _b 


0.52, 1.92 


1.9 x W~ b 


0.54, 1.87 


12.42 


Float 


2 x 10~ 5 


0.52, 1.92 


2.0 x 10^ 


0.54, 1.87 


12.06 


Fix 


3 x 10~ 5 


0.52, 1.92 


3.0 x 10^ 


0.50, 2.01 


12.17 


Float 


3 x lO" 5 


0.52, 1.92 


3.0 x 10~ b 


0.51, 1.94 


11.88 


Fix 


7 x 10~ 5 


0.52, 1.92 


7.2 x W~ b 


0.45, 2.23 


10.73 


Float 


7 x 10~ 5 


0.52, 1.92 


7.0 x 10~ b 


0.61, 1.63 


9.74 


Fix 


1 x 10- 4 


0.52, 1.92 


1.0 x 10" 4 


0.51, 1.94 


12.96 
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1 x 10" 4 


0.52, 1.92 


1.0 x 10" 4 


0.58, 1.71 


11.71 
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FIG. 7. The contours that arise from including the flux 
uncertainty in the x 2 . 



a few brief comments are possible. First of all, in the 
preceding analysis, we have neglected the contribution 
from the backgrounds. In the reactor experiment, there 
is a distinctive delayed neutron-capture signature, which 
results in a substantial reduction of backgrounds. The 
background level was conservatively estimated to be 
20:1 j|] (i.e., probably an overestimate). It is not so 
unreasonable to ignore the backgrounds. Moreover, one 
expects that backgrounds will be a relatively steeply 
falling function of energy. So, if the energy where the 
interesting oscillation effects occur is high enough, then 
the assumption of no backgrounds is even safer. This 
happens at the larger Am 2 values, as can clearly be seen 
from Eqn. (2). In addition, KamLAND hopes to get 
a handle on the backgrounds by utilizing the fact that 
power usage, and hence neutrino flux, varies seasonally 
H . The amount of power produced over time will be 
available from the power companies. 

We would like to briefly mention other possible sources 
of systematic uncertainties. As mentioned earlier in 
the paper, there is an issue of the energy resolution of 
the detector. This could smear out the location of the 
"dip" that is so nicely seen, in Fig. ||(a). However, the 
resolution expected is better than the size of the bins 
and hence is not expected to affect the results. A more 
difficult issue would be to accurately callibrate the energy 
measurement of the detector of this size. There are plans 
to do so; see ||. There will be a source of backgrounds 
from geological neutrinos in the 2-3 MeV region which 
may or may not be important. However, a simple energy 
cut should remove them completely. 



VIII. CONCLUSIONS 



The KamLAND reactor neutrino experiment should 
allow for an accurate measurement of Am 2 , as long 
as the solution to the solar neutrino problem is in 



the LMA region. Incomplete knowledge of the fuel 
composition at the reactors will represent an important 
source of experimental error. In this letter, we assumed 
perfect detector efficiency, and that backgrounds were 
well understood. It is interesting to note that the 
precision to which KamLAND is able to measure Am 2 
may effect the ability of a muon source neutrino factory 
to extract information about the CP violating phase in 
the MNS matrix. CP violation is an phenomenon that 
requires at least three generations, the effects of the sub- 
leading oscillation is crucial, and the matter effect needs 
to be accurately subtracted in order to determine Sep- 
The measurements at KamLAND should prove crucial 
for this purpose. 
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NOTE ADDED 

During the final preparation of this manuscript, we 
became aware of papers on a similar topic p3| . [l4| . 
Both of them discussed how accurately KamLAND will 
determine the oscillation parameters. Ref. |l3[ ] fur- 
ther included the three-generation mixing effects, while 
Ref. g studied Am 2 > 2 x 1(T 4 eV 2 as well. Neither 
of them discusses the effect of fuel composition and 
overall flux normalization quantitatively, however. Other 
conclusions appear to be consistent. 
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